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According to quantum field theory, the inflaton equation of motion does not have the local form that is gen- 
erally assumed for cosmological purposes. In particular, earlier investigations of the nonequilibrium dynamics 
of an inflaton that decays into scalar particles suggest that the loss of inflaton energy is not well approximated 
by the local friction term derived from linear response theory. We extend this analysis to the case of an inflaton 
that decays into fermions, and reach broadly the same conclusion. 
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I. INTRODUCTION 

In the study of inflationary cosmology, it is often assumed 
that a classical inflaton field <j> obeys a local equation of mo- 
tion of the form 



+ 7y(0)0 + y'(0) =0, 



(1) 



where the friction term ri{(j))(f) accounts for the loss of in- 
flaton energy by radiation of particles. (We omit the usual 
Hubble damping term SHcf), for reasons discussed below.) 
In particular, "warm inflation" scenarios (see, for example, 
||Jj,2b^^ 5]) envisage a friction coefficient r]{(f>) large enough 
for the inflaton motion to be perpetually overdamped, thus 
avoiding the usual process of reheating. 

However, an equation of this form does not arise automat- 
ically from quantum field theory. Indeed, if the inflaton cou- 
ples biquadratically to scalar fields x, and via Yukawa interac- 
tions to spin-i fields tp, its equation of motion is generically 
of the form 



^ + r(</))+.gi(x')0 + 52(^^) 



= 



(2) 



where the expectation values (x^) and are nonlocal 

functionals of (p. That is, they depend on the entire history 
of <j), at all times prior to the time t at which the expectation 
values are evaluated. The issue we address here is whether 
these expectation values can be adequately approximated by 
a local friction term, together with a local correction to the 
effective potential V{(j>), leading to the equation of motion 
([T]i. To focus on this issue, which is substantially independent 
of the background metric, we simplify matters throughout by 
studying quantum field theory in Minkowski spacetime, with 
H = 0. 

If the motion of <j) is slow enough, it seems intuitively rea- 
sonable that suitable approximations to these expectation val- 
ues can be obtained by considering the linear response of a 
state of thermal equilibrium to a time-dependent perturbation, 
4>{t). Several implementations of this general idea can be 
found in the literature (see, for example lflll6l[7i[8ll9l [Toll ). and 
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lead to well-defined estimates of the friction coefficient 77 (</)). 
However, the assumption that this problem can adequately be 
treated on the basis of linear perturbations about equilibrium 
thermal field theory deserves further investigation. If this as- 
sumption is valid, then the results should be reproduced in the 
slow evolution limit of nonequilibrium field theory. 

For the case of an inflaton coupled to a scalar field, this does 
not seem to be true ifTTl 1 1311 . Formally, one can derive 
approximate, local nonequilibrium evolution equations, but a 
time-derivative expansion of the solution to these equations, 
which is needed to arrive at ([T]i, does not exist. With the fur- 
ther approximation of replacing non-equilibrium self-energies 
with their equilibrium counterparts, a time-derivative expan- 
sion becomes possible, leading to what we call the adiabatic 
approximation, which essentially reproduces the results of lin- 
ear response theory. However, numerical investigation shows 
that the motion generated by the nonequilibrium equations is 
not well approximated by the adiabatic approximation. 

These results might well be specific to the particular field 
theory considered. In this paper, we investigate this to some 
extent by extending the analysis of ITll [T2I u3\ to the case of 
an inflaton that decays into spin-i fermions. In principle, the 
calculations are quite similar, but the technicalities associated 
with spin-i fields are slightly more complicated. 

In section we derive exact evolution equations for the 
number densities Nkit), and auxiliary functions Vk{t), which 
together specify the nonequilibrium state of fermions that are 
free, except for a i/i-dependent mass arising from a Yukawa 
coupling to the classical inflaton. These are the fermionic 
analogues of the exact equations derived for scalar fields by 
Morikawa and Sasaki [7]. As in the scalar-field case, they do 
not describe any significant frictional effect, because particles 
radiated by the inflaton motion can easily be reabsorbed. To 
achieve a permanent transfer of energy from the inflaton to 
the system of particles, these particles must interact with each 
other, and dissipate energy by scattering. 

An approximate set of nonequilibrium evolution equations 
is obtained in section |llll by adapting the analysis described 
by Lawrie and McKernan in [14]. The closed-time-path meth- 
ods used to develop this approximation are very different from 
the straightforward (if long-winded) operator methods avail- 
able for free field theories, and the fact that the equations 
of section are recovered in the free-field limit is a crucial 
check on this method of approximation. We again find that 
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the nonequilibrium evolution equations do not admit a solu- 
tion in the form of a time-derivative expansion (in fact, this 
seems to be a general, model-independent result), but an adi- 
abatic approximation can be devised, as explained in section 

Numerical solutions to the nonequilibrium evolution equa- 
tions, and to the adiabatic approximation to these equations, 
are described in section |IV] As in the previous case of an in- 
flaton that decays into scalars, we find that these solutions do 
not agree well. In general, the adiabatic approximation ap- 
pears to overestimate the damping effect of particle radiation 
on the inflaton motion; in particular, we find that overdamping 
predicted by the adiabatic approximation, is not exhibited by 
the nonequilibrium evolution. 

A self-contained summary and discussion of these results is 
given in sectionlV] 



n. FREE-FIELD EVOLUTION EQUATIONS 

In this section, we derive a set of exact equations for the 
time evolution of a system of fermions that are free, apart from 
a coupling to the classical inflaton; this gives rise to an effec- 
tive time-dependent mass m{t) = + g4'{t), where 
denotes the bare mass. These equations serve as a valuable 
check on the validity and interpretation of the approximate 
evolution equations we obtain in section |III] for the system of 
interacting fermions. 

We look for a formal solution of the Dirac equation with 
time dependent mass. 



with 



[z7^5^ - m{t)] ^{x) = 



(3) 



by expanding the field ijj{x) = ip{x,t) in terms of two 
positive-energy spinors u± and their charge conjugates — 



(2^ 



[Ak (t; i)u+ {k, i) + Bk [t- i)u^ (fe, i) 



+Ckit; i)u%{-k, t) + Dk{t; i)u<L{-k, t)]. 



(4) 



The spinors u± (fe, i) are helicity eigenstates, with helicity ^h, 
h — ±1, as indicated by the subscript and, at an arbitrary ref- 
erence time i, are solutions of the instantaneous Dkac equa- 
tion 



+j-k + m u±{k, i) = 0, 



(5) 



where m — m{i) and Clk = y/k^ + m^, normalized such 
that ul^{k,i)uh'{k,i) = Shh'- Dependence on the time t 
is contained in the functions Ak{t; t), Bk{t; t), . . . and these 
coefficients also depend on i, in such a way that ip{x) itself 
is independent of t. By substituting the Ansatz dU into the 
Dirac equation (|3]l, we find equations for the t dependence of 
Ak{t\ i), Bk{t; i), . . .. The change of variables 



Ak = CkAk - SkCk Bk = CkBk - SkDk 
Ck = SkAk + CkCk Dk = -SkBk - CkDk 



(6) 



Ck 



flk + fn 



\k\ 



Sk 



2{lk{Clk + in) 



(7) 



(satisfying cj. + = 1) brings these equations into a standard 



form used in 111411 . namely 

idt 
idt 

The matrix 



Ak{t;i) 
Ck(t;i) 

Bk{t-i) 

Dk{t-i) 



Tk{t) 
Tkit) 



Ak{t-i) 
Ck{t;i) 

(Bk{t\i) 
\Dk{t;i) 



(8) 



m{t) — fc 



(with k 



Tkit) = 

|fe|) has eigenvalues ±17fc(i), where ilk{t) 



(9) 



•\/ m{t)'^ + fc^. Correspondingly, two orthogonal solutions to 
the equation idtF{t) = Tk{t)F{t) can be written as 



.(+) 



{t;i) 



= (fk{t\i)\^rj-^^tJlT^(t')dt' 



Ck 
-Sk 



(10) 



v9fe(^;t) 

where T denotes time ordering. At the reference time i, pj^^ 
and f[ ^ are positive- and negative-energy solutions respec- 



(±) 



tively (that is, idtF^ 

The general solution for Ak{t: i), Bk{t; i), . . . can thus be 
written as 



B{t;i) 
D{t;i) 



b^{k,i) 



9k{t\t) 

.fk{t;i) 
9k{t;i) 



d\{-k,t) 



-9*k{t;i) 
fkit-^i) 



(12) 



and the canonical equal-time anticommutator of ij^ix) and 
■4>\x) implies that the coefficients 6/i(fe, i) and dh{k, i) have 
the usual anticommutator algebra 

{bh{k,t),bl{k',i)} - (2^)35(fc-fc')<5w,etc 

for particle creation and annihilation operators. 

To find how these operators depend on t we make use of the 
fact that dfTpix) = 0. On substituting Q and ( fT2] i into (|4|, we 
find (after some algebra) that 



dMk, i) = -ihkhhik, i) + hXkdli-k, t) 
didh{k,i) = -iClkdhikJ) - h\kbl{^k,t) 



where 



A. = 



kdfQk 



kd^rh 



(13) 



(14) 
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In a translation- and rotation-invariant state, expectation val- 
ues of bilinear products of creation and annihilation oper- 
ators are characterized by functions Nh{k,t), Nfi{k,t) and 
iyh{k, t), defined by 

{bl{k,i)bh'{k',i)) = {2TT)H{k--k')5hh'Nh{k,i) 

{di{k,i)dh'{k',i)) = {27TfS{k - k')6hh'Nh{kJ) (15) 
{dhi~k,i)bh'{k',i)) - {27rf S{k - k')Shh'Mk,i). 

Because the mode functions -F';^^^ {t, i) are positive-frequency 
only for times t near t, we can interpret Nfi{k,t) as the 
number density of particles present at time t with momen- 
tum k and helicity h, and Nfi{k, t) as the corresponding den- 
sity of antiparticles. The functions Vh{k,t) measure the off- 
diagonality of the density matrix in the basis specified by the 
modes F^^\ and have no classical analogue. 

The expectation value (i/'(a;, t)ip{x, t)) is conveniently ex- 
pressed in terms of the number densities at time t, and we 
achieve this by setting t = t. We will restrict attention to states 
in which N+{k,t), N-{k,t), N+{k,t) and N-{k,t) are all 
equal, say to Nk{t), while l'-^-{k,t) = —v^{k^t) = Vk{t)- It 
is easily checked using ( fT3T l that these equalities are preserved 
by the time evolution, and that the real number density Nk [t) 
and the complex function Vk {t) obey the evolution equations 



by Lawrie and McKernan in fT?]. We follow their derivation 
closely, except that we introduce, as in sectionHH an evolving 
reference time t, which allows us to make explicit the effect 
of particle creation. 

Within the closed-time-path (CTP) formalism (see, e.g. 
lITsi [T6ll and references therein), we seek an approximation 
to the full two-point function 

{ {M,x)i^p{t' ,x')) {f[,p^{t,x)i>p{t',x'))l 

(19) 

from which the expectation value ('(/'(a;, t)tlj{x, t)) can be ob- 
tained as a special case. (Here, a, f3 are Dirac indices, while 
a,b = 1,2 distinguish the indicated operator orderings.) It is 
shown in (l^ that this can be expressed (after a spatial Fourier 
transform) in terms of a single 4x4 matrix 7i(t, t'; k) as 



5("'')(tt';fe) = 



'H{t,t';k) ■H{t,t';k) 

yH{t,t';k) ■H{t,t';k) 

^■H{t,t';k) 'H{t,t';k) 

(Hit.t'-k) n{t,t'-k) 




dtNkit) 



krn{t) 



fit) 



(16) 



(17) 

where is the real part of i^. The expectation value of inter- 
est is given in terms of these functions by 

d^k 



For any matrix M, we define M = '^°M^-f^ and M = 
[C~^7°A/7"C]"^ where C is the charge conjugation matrix 
and ^ indicates the transpose. 

As explained in detail in ifTil] (and in lfT2l [TtIi for the case 
of (/) coupling to a scalar field), we construct a partially- 
resummed perturbation theory as follows. Given a CTP ac- 
tion / = + where /'^^ is quadratic in the fields and 
/(^^^ is of higher order, we take its lowest-order part to be 



(V'(a;,t)^(a;,i)) =4 



{2nf2^k 



[m{2Nk - 1) + 2kv^] , /o(V>) = + / d^xi^aMabi^b = / d^X^PaVab^b 



(18) 

with all the quantities to, ilfc, Nk, Vk evaluated at time t. The 
evolution equations ( fTSI l and ( [TtT i are the fermionic analogues 
of those derived by Morikawa and Sasaki iQ] for the case of 
the inflaton field coupled to a scalar field. The right-hand side 
of (fTSI l is a particle 'creation' term, proportional to d(p/dt, 
but it can be either positive or negative as ly^ oscillates, and 
a numerical solution shows that there is no permanent trans- 
fer of energy from the inflaton to the system of fermions. To 
see a frictional effect, it is necessary to introduce interactions 
between the fermions so that, roughly speaking, the created 
particles lose coherence with the inflaton, and a permanent 
dissipation of energy results. 



III. EVOLUTION EQUATIONS FOR AN INTERACTING 
SYSTEM 



A. Derivation 

An approximation scheme for studying the non-equilibrium 
evolution of a system of interacting fermions was developed 



(21) 

where A4 is a counterterm to be determined, and treat the re- 
mainder, /int = /(>2) _ J d'^x'ipaMab'^b as a perturbation. 
(For a = 1,2, tpa are the independent path-integration vari- 
ables that inhabit the real-time branches of the closed time 
path; we have indicated explicitly only the contribution of one 
spin-i field.) The lowest order propagator S{t,t') (we sup- 
press the momentum argument k) obeys the equation 

= lSabS{t-t') (22) 

where the differential operator T), defined by ( |2TI) . is 

I5(t,a.)=(^°(f*) _^^l^g^^)+Mit), (23) 

and Vq = i^^dt — 7 • fc — m{t) is the usual Dirac operator. 
Given suitable constraints on the form of the counterterm M., 
the solution for S'^'^''^ {t, t') can be written in the form (|20] i in 
terms of a function H{t, t') which is the lowest-order approx- 
imation to 7Y(i, t'). 
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In this way, we obtain a reorganized perturbation theory, in 
which the vertices impHed by /(>2) are supplemented by the 
counterterm . In particular, self-energies have the form 



S„fc(tt') = Mab{t)S{t - t') + Y^'Tit^t') 



(24) 



and we can optimize the propagators 5^"^^ as an approxima- 
tion to the full 2-point functions 5'-''''^ by choosing Mab to 
cancel a local part of the loop diagrams in ■ 

The solution of ( |22] | is explained in detail in ifTill . The ma- 
trix H[t^ t') from which S^°-^^ (i, t') is built can be expressed 

in terms of the 7 matrices, the mode functions F^'^ (t, i) in- 
troduced in section HH and two auxiliary functions Ni^{t,t) 
and Afc {t, t) which, in the free-field limit, reduce to Nk (t) 
and ~Vk{t) respectively. These functions satisfy differential 
equations in the variable t, which arise from ( |22] |. In ifTill . the 
reference time i was taken to be fixed at at the instant when 
an initial state is specified, but here we find it convenient to 
allow t to vary, and eventually, as in section |II] to set t = t. 
The dependence of Nk{t,t) and Ak{t,t) on t can be found 
from the fact that 5*^°^-' {t, t') must be independent of i. 
We now define 

Nk{t) = Nk{t,t) and Vk{t) = -l^k{t,t). (25) 
Evidently, we have 

dtNu (t) = [dtNk {t, i) + d-.Nk [t, £)] (26) 

and similarly for Vkit)- Applying the results of [Ti!], we find 
that these two equations have the form 

dtNk{t) = -^Tk{t)[2Nu{t)-l]+ak{t)+^^^v^{t) 



krn{t) 



[2Nk{t) - 1] . 



(27) 



(28) 



The functions Tk (t), which can be interpreted as a quasiparti- 
cle decay width, and ak {t) appear in the counterterm A4 and, 
according to the strategy outlined above, are to be determined 
from the self-energies 

A concrete realization of this strategy is described in 
for the case that the fermions interact through a Yukawa cou- 
pUng to a scalar field x(a;, i) of mass M, 



(29) 



and we consider the same model here, assuming for simplic- 
ity that X does not couple to the inflaton, so M is constant. 
Evaluating the relevant self-energies to 1-loop order gives 



ak{t) 



647r2 

X [rip + Nk' 



(30) 



647r2 ' J ilkilk'tup 

x[npil-Nk')~il + np)Nk']. (31) 



Here, m stands for the fermion mass m{t), np{t) is 
(approximately-see below) the number density of scalars, 
p = fc + fc' is the momentum of a scalar produced, say, 
by the collision of two fermions of momenta k and fc', and 
LUp — + A/2. In this approximation, the first two terms 

in equation (|27] | combine to give 



stit) 



-Tk{t) [2Nk{t) - I] + ak{t) 



i2lT-^ J ilkilk'^p 



X [np{l^Nk){l-Nk')-il + np)NkNk'] (32) 

which we recognize as the Boltzmann scattering integral cor- 
responding to the decay of a scalar into two fermions and 
the inverse production process. The evolution of np{t) is de- 
scribed by a similar Boltzmann-like equation 



dtn. 



Qk 



Qk') 



IGtt^ '' J Qk^k'^p 

X [np{l-Nk)il-Nk')-il + np)NkNk'], (33) 

where now k ^ p—k'. This contains no particle creation term 
(analogous to the last term of ( l27b ) owing to our assumption 
that x{x,t) does not couple to the inflaton. 

At this point, we have a closed system of evolution equa- 
tions, consisting of dzTl ). ( |28] ). ( |33] ) and the inflaton equation 
of motion, which we now take to have the form 



tru 



= 0. 



(34) 



Here, we have approximated {■ip{x,t)ip{x,t)) by its lowest- 
order contribution, — (27r)^^/ d^kTTH{t,t' ,k). It consists 
of a local function of <j){t), which contributes to the effective 
potential, together with the term 



(l/"/')trunc = 4 



d^k 



{2^f^k{t) 



[m{t)Nk{t) + kv^{t)\ 



(35) 

Since we are concerned with the frictional effects arising from 
('0V')trunc, we simplify matters by replacing the whole effec- 
tive potential with ^m^^^. For later use, we observe that the 
total energy given by 



dkk' 
~2^ 



[iOkUk + AVLkNk] (36) 



is exactly conserved by this system of evolution equations. 

In the interest of accuracy, two remarks are in order First, 
the representation ( l20l i of the full 2-point functions is valid 
for a CP-invariant theory. In practice, this means only that 
CP-breaking interactions will not be resummed by the coun- 
terterm Ai. Second, the quantities denoted in this section by 
Nk{t), Vkit) and np{t) are functions that arise in the process 
of solving the appropriate differential equations for spinor and 
scalar propagators. Only in the free-field limit can they be un- 
ambiguously identified in terms of expectation values of cre- 
ation and annihilation operators as in (flSl l. 
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B. Renormalization 



C. Adiabatic Approximation 



The integral dSSl) that appears in the non-equilibrium equa- 
tion of motion ( [34] i has a divergent contribution from the term 
i/^ in the integrand. To isolate this divergence, consider the 
evolution equation ( |28] ) in the limit that k is large. Assuming 
that the thermal terms proportional to and can be ne- 
glected (which is to be expected by analogy with equilibrium 
field theory, see e.g. section 3.5 of Ref II18I1 . but hard to prove 
here) we find 



With the approximations developed above, the functions 
Nk{t) and i^fc(t) that appear in {ipip) trunc obey local evolution 
equations (|27| i and ( |28] l, but the solution to these equations is 
non-local, i.e. (V'V')tiunc evaluated at time t depends on (/)(t') 
at all times t' prior to t. One can attempt to derive a local 
approximation of the form 



5('0V')tri 



Ay'(0) + 7/(0)0, 



(43) 



The solution can be written as 

/i\ —2ikt 

Vk{t) = e 

and on repeatedly integrating by parts we find 



(37) 



(38) 



—2ikt 



-(O) + ^^(O)-^0(O) 



(39) 



Now, contributions to v^it) that vanish no faster than k^^ as 
k ^ oo will yield a divergent integral in ( [35] ). To avoid this, 
we firstly choose initial conditions which have the form 



Mo) = -^m + ^m + o{k-') 



8fc3 



(40) 



for large k, so that the square bracket in ( |39l ) vanishes. We are 
then left with the term {g/8k^)(j){t) which can be canceled 
by a wavefunction renormalization (p Z^^'^tj) in the equa- 
tion of motion ( l34l i. This is, in fact, the same wavefunction 
renormalization that would be required in a zero-temperature 
quantum field theory with the Yukawa coupling ( |29] l. With a 
suitable choice for Z, the renormalized integral is 



^ 4 



2 r 



dkk 
gk 



m{t)Nk{t) + kiyj^jt) 



8(P 



(41) 



where mp is a constant. The corresponding renormalized ex- 
pression for the conserved total energy is 



2 i2 



dkk^ 



i^knk + 4:ilkNk 



16(fc2 



(42) 



Note that, since vj^ has a contribution of order k^^, the spec- 
trum of particles created by the last term in dZTl l has a cor- 
responding contribution of order fc""*, and so therefore does 
Nk- This gives a convergent integral in i4T[ . but a divergent 
contribution to £''°", which is cancelled by the last term in 
dUli. 



which depends only on the values of and at time t, by 
means of a time-derivative expansion, but we show below that, 
as in the case of an inflaton that decays into scalar particles 
flT], this does not work. 

To facilitate a time-derivative expansion, we introduce a 
formal expansion parameter e multiplying time derivatives 



edtNkit) 



St 



edtvkit) = -2[ink{t)+Tk{t)H{t) 
gkdt(j){t) 



2nl{t) 



[2Nk{t) - 1] 



(44) 
(45) 

(46) 



where we have also used m{t) — gdt4>{t). We now expand 
Nk{t), np{t) and i^fc(t) in powers of e around the equilibrium 
distributions iV^''(t), n^'i(t) and j^l^{t) 

Nk{t) = Nl\t)+e6Nk{t) + Oie^) (47) 
np{t) = n'pit) + eSupit) + 0{e^) (48) 
^fcW = 'yr{t) + edi^k{t) + 0{e^). (49) 

Substituting this expansion into equations (|44|)-(|46]| we find to 
leading order that i^'j,'^ = 0, while 



1 



(50) 



for some inverse temperature /3 are the usual Fermi-Dirac and 
Bose-Einstein distributions that make the Boltzmann scatter- 
ing integrals vanish. At order e, SNk and Sup are solutions of 
the integral equations 



dtN? 



dk' ^^'^ 



= / dk' 



SNk' 



SNk' 



SNk' 



dp' " 



dp' 



Supi 
SS^ 



5n^ 



Sup' 



ip' 



(51) 



where |cq means that we set N = and n — n'"^ after dif- 
ferentiation. However, these equations are not self-consistent 
and therefore have no solution. It is easily seen from ( l32b and 
(O that 

S = 2 j d^kSt + j d^pS^ = 0. (52) 
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This reflects the fact that if TV = J (Pk Nk is the total number 
density of fermions of each hehcity (and also the number of 
antifermions of each helicity) and n — J d^p Up is the total 
number density of scalars, then the quantity 2N + n is con- 
served by the processes x ^ V'V'- (For example, one fermion 
of a given helicity is produced, on average, in every two x de- 
cays.) In the presence of particle production, this conservation 
law no longer holds, but the scattering integrals still obey the 
sum rule {5% . The integral equations ( BTI ) therefore entail 



2dt J d^kN^+dt J d'pnl'i 



= 0, 



dp oririi 



(53) 



but this is not true for the equilibrium distributions ( |50l l. Thus, 
although the evolution equations themselves are consistent 
(albeit approximate), they do not admit a solution in the form 
of a time-derivative expansion. In principle, therefore, we 
cannot approximate (V'V')trunc in the form ( |43] i or the infla- 
ton equation of motion in the local form ([T]) by means of a 
time-derivative expansion. 

A local approximation to the equation of motion can how- 
ever be obtained by resorting to a further approximation, 
which consists in replacing Nk{t) and np{t) with their equi- 
librium values (ISOl l for the purpose of evaluating the functions 
Ffe, which we then denote by T^^, and ak- When this is done, 
the 0(e) equations for 5Nk and dvk are 



(54) 



= -2[zl7, + Tl'']Si^k - X [2^r - 1] -(55) 

On rearranging these, we find that the inflaton equation of mo- 
tion can be written as 



(56) 



with 



dkP 
dkk^ 



m{t)Nl\t) + kv^'-^\t) 
m(t) 



and 



m{t)SNk{t) + kSiy^{t)] 



(57) 



dkk 



9 
47r2 



dkk^ 

cc 

fe 



[1 - 2iV°^(i)] }(58) 



In the zero temperature hmit (N^'^ = 0) and taking << fl^ 
this expression agrees with that found in [lO] (apart from a 
factor of 2 which we have not been able to account for). We 



refer to the approximation represented by (I56]l-(l58ll. together 
with a suitable prescription for determining the temperature 
(3^^, as the adiabatic approximation. 

To estimate time evolution over any extended period of 
time within the adiabatic approximation, we must allow (3 to 
change in an appropriate manner We obtain a suitable pre- 
scription by requiring that the energy (|36] l be conserved by 
the adiabatic evolution when Uk and Nk have their equilib- 
rium values. 



IV. NUMERICAL INVESTIGATION 

The approximate nonequilibrium evolution equations de- 
veloped in section IIII A I lead to the local equation of motion 
(|56] | only when we make the extra approximation of replac- 
ing self-energies (from which the functions ak and Tk are 
derived) with their values in a state of thermal equilibrium, 
together with the time-derivative expansion which then be- 
comes possible. On the face of it, one might expect these fur- 
ther approximations to be fairly harmless, at least for a system 
that does not evolve too fast, and does not stray too far from 
thermal equilibrium. However, the numerical calculations de- 
scribed in 1 13] for the case of an inflaton coupled to scalar 
particles showed that the nonequilibrium evolution governed 
by a set of equations analogous to those of section IIII A I is 
not well represented by the corresponding adiabatic approx- 
imation, even under circumstances when that approximation 
might seem to be reasonably well justified. 

We therefore discuss in this section a similar numerical 
comparison of the time evolution generated by the nonequilib- 
rium equations described in sections UlI Al and [ III Bl with that 
generated by the adiabatic approximation of section UlI CI To 
discretize the nonequilibrium time evolution, we use a semi- 
implicit method defined schematically by 



dx 




dx 




dt 


n 




n+l 



5t 
2 ' 



(59) 



where x represents the whole collection of variables 0, Nk, 
Uk, and Vk and n labels time steps. In the absence of the 'ipipx 
interaction, it is straightforward to show that this method is 
stable for our problem. Crucially, it is also simple enough that 
we can implement a discretized version of the renormaliza- 
tion scheme of section IIII B I that is exact up to rounding er- 
rors. When interactions are included, stability of the resulting 
set of integro-differential equations is very hard to analyse. 
We do in fact see evidence of instability when the calcula- 
tions are pursued for sufficiently long times, but the results 
presented below are for periods of time for which the effects 
of this instability appear to be negligible. Moreover, although 
the renormalized total energy (l42l) is exactly conserved in the 
continuum limit, there is no local discretized version of this 
expression that is exactly conserved by the discretized evolu- 
tion equations. As a rough check on our computations, we 
have evaluated a discretized version of dill, but we do not 
expect it to be exactly conserved. In particular, divergent con- 
tributions to the integral, which cancel by construction in the 
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continuum limit, are not guaranteed to cancel when time evo- 
lution is discretized. In practice, all the initial conditions we 
have considered lead to an oscillatory solution for (j){t); we 
find that energy acquires an oscillatory component (with am- 
plitude of order 5% of the total), consistent with incomplete 
cancellation of the divergent integral proportional to 0^ , to- 
gether with a drift amounting to a few percent over the time 
intervals of interest. 

The adiabatic approximation is not affected by these issues. 
Here, no renormalization is required; semi-implicit discretiza- 
tion of the one remaining differential equation, the equation of 
motion ( l56b . leads to a stable problem; and exact conservation 
of energy is imposed as the criterion for determining f3. In- 
tuitively, this approximation ought to be good when is 
slowly varying and thermalization of the created particles is 
efficient. A rough criterion for efficient thermahzation is 



or. 



T I dm / dt\ <C m 



(7 = g4>/mTk=a < 1 



(60) 



(61) 



where r is a relaxation time, which we take to be r = 1 /Tk=Q- 
To obtain slow variation in we would ideally like to 
identify a set of parameters for which the motion is over- 
damped. Within the adiabatic approximation, one way of 
achieving this is to make the first term in the friction coeffi- 
cient (ISSl l large, by making Tk relatively small, but this tends 
to make thermalization slow, at least according to the crite- 
rion f6i} . The parameter space of couplings, masses and ini- 
tial conditions is, even for the simple model considered here, 
too large to permit a systematic exploration, but we have not 
succeeded in locating a set of parameter values that leads both 
to overdamped motion and to a very small value of a. (This 
was also true of the investigation reported in for an infla- 
ton decaying into scalars; in the context of warm inflation, it 
seems to be generally hard to devise simple field theory mod- 
els that exhibit efficient thermalization as well as other desir- 
able properties |4, 9, 19].) 

Nevertheless, we can find a situation in which the adiabatic 
approximation predicts overdamped motion, and a is not too 
large. Figure [T] shows the evolution of (f){t) calculated from 
the nonequilibrium evolution equations and from the adia- 
batic approximation, with the parameter values g = = 1- 
TO^ — ra^ — 3, m^p — 0.1, and initial conditions 0(0) — 1, 
(/)(0) = 0, /3 = 1. The solid curve is the motion gener- 
ated by the adiabatic approximation, which is indeed over- 
damped; we see, however, that the motion generated by the 
non-equilibrium equations without the adiabatic approxima- 
tion (dashed curve) although it exhibits some damping is far 
from overdamped. 

As in [Ts], this striking discrepancy does not seem to be ac- 
counted for by inefficient thermalization. In figure |2] we show 
the early stages of the motion, along with the parameter a as 
estimated within the adiabatic approximation. During this pe- 
riod of time, a remains below 1. Thus, within the adiabatic 
approximation, one might well conclude that this approxima- 
tion is reasonably self-consistent. While thermalization is not 




FIG. 1 : Evolution of in both the adiabatic and non-equilibrium so- 
lutions. The initial conditions predict overdamping for the adiabatic 
solution. 
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FIG. 2: Evolution of (f> and a for early times. The thermalization is 
relatively efficient in this period and the two solutions are still very 
different 



efficient enough to guarantee good accuracy, one would prob- 
ably not anticipate the gross discrepancy between the two ap- 
proximations for (t){t). Moreover, we compare in figure[3]the 
function Ffe calculated at t « 3 from the nonequilibrium num- 
ber densities with that calculated in the adiabatic approxima- 
tion from an exactly thermal distribution. While some dif- 
ference is visible, it seems clear that the small quantitative 
difference cannot in itself account for difference between the 
overdamped and oscillatory motions apparent in figure|2] 

In fact, the difference between these motions can be under- 
stood in the way suggested in ifTsll . At < = 0, we chose an 
initial state of exact thermal equilibrium. Denote by A{(f)) 
the integral ( [35] ) evaluated in this state, and let U{(l)) = 
+ A{(f)). At early times, when the value of has not 
changed by much, the nonequilibrium equation of motion ( |34| ) 
is given approximately by 



+ C/(0)«O, 



(62) 
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FIG. 3: Momentum distribution of Ffe at f « 3 for both the adiabatic 
and the non-equilibrium solutions. 



while the adiabatic approximation to this equation is 

^ + U{(t>) + w 0. (63) 

As seen from figure|2l where a cx \(t>\, the velocity cj) in the adi- 
abatic approximation quickly approaches the "terminal veloc- 
ity" K, —U {(/))/ ri{(p) and thereafter </> and evolve slowly, 
with (j) Ri 0. According to the original equation, by contrast, (p 
is still approximately equal to —U{4)), leading to the dashed 
curve. Clearly, in these circumstances, the original equation of 
motion is not well represented by the adiabatic approximation, 
and this is substantially independent of how well the number 
densities used to calculate U {(p) and 77(0) are thermalized. 

An example of the more generic case of underdamped mo- 
tion is shown in figure |4] We see that the nonequilibrium evo- 
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FIG. 4: Underdamped motion of (j> for both the adiabatic and non- 
equilibrium solutions 

lution shows significant damping, but that the rate of damping 
is overestimated by the adiabatic approximation. This differs 
somewhat from the case of inflaton decay into scalars [IS]. In 
that case, the adiabatic approximation leads to a friction coef- 



ficient roughly proportional to cjp, and thus to an overestimate 
of damping at large amplitudes, but an underestimate at small 
amplitudes. 



V. SUMMARY AND DISCUSSION 

In this paper, we have sought to address a very specific is- 
sue: can the expectation value {ij^ip) appearing in the generic 
inflaton equation of motion dU, which is a non-local func- 
tional of (j), be adequately approximated by a local function 
of (f){t) and 4i{t), leading to a local equation of motion such 
as ([TJ. The conclusion suggested by the results of the two 
preceding sections is that it cannot be, but these results are 
neither exact nor wholly rigorous, so it is worthwhile to sum- 
marize carefully the line of argument we have adopted. 

Derivations of the local friction coefficient 77(0) to be found 
in the literature vary somewhat in the details of implementa- 
tion, but generally involve the assumption, characteristic of 
linear response theory, that under suitable conditions the state 
of the particles can be treated according to equilibrium field 
theory. The possibility of obtaining a local equation of mo- 
tion then depends on the behaviour of the relevant integral 
kernels, which can be investigated in perturbation theory for 
specific models (see, for example [lOJ and, for recent discus- 
sions ll20ll2lll '). We seek to avoid this assumption by studying 
instead the slow-evolution limit of the nonequilibrium state. 
In principle, this is certainly the more correct procedure; in 
practice, however, nonequilibrium evolution is extremely hard 
to deal with by analytic means. 

The approximation scheme sketched in sections IIII Al and 
IIIIBI is the lowest non-trivial order of a partially resummed 
perturbation theory. It leads to a local set of evolution equa- 
tions, but we find that these equations do not admit a solution 
in the form of a time-derivative expansion, which is needed to 
derive the local equation of motion ([1]), and consequently that 
this local equation of motion does not, in principle, exist. We 
cannot altogether rule out the possibility that this result is a 
feature of our approximation scheme, rather than of the exact 
dynamics. Quantitatively, indeed, the accuracy of the approx- 
imation is likely to be quite modest, though we have no good 
way of assessing this. However, the fact that it reproduces the 
exact dynamics of section|lI]in the free-field limit, and that the 
Boltzmann-like scattering terms have a sensible kinetic-theory 
interpretation, offers some reassurance that the approximation 
captures essential features of the nonequilibrium dynamics. 
Moreover, any improved method of studying nonequilibrium 
dynamics (for example, the 2-particle-irreducible formalism 
developed by Berges and coworkers; see 1^25] and references 
therein, especially [23]) is intrinsically non-local. It seems 
very likely that any attempt to derive a local equation of mo- 
tion must follow a sequence of approximations more or less 
equivalent to that used here, and would fail at the same point. 

Even if the local equation of motion is formally not valid, 
it might nevertheless be a good approximation under suit- 
able circumstances. Indeed, we can recover essentially the 
same result for !]{([)) as that given by linear response methods 
by taking the further step of replacing nonequilibrium self- 



9 



energies with those calculated in a state of thermal equilib- 
rium: we called this the adiabatic approximation to our non- 
equilibrium evolution equations. Intuitively, this would seem 
to be a fairly innocuous step if the evolution is fairly slow, 
and thermalization is not too inefficient. To test this, we ob- 
tained numerically, in section IIVI solutions for the motion of 
(f) both with and without the adiabatic approximation. In order 
to focus on the frictional effect (and to simplify the numerical 
calculation somewhat) we took the whole effective potential 
(including the local part of {"ijjii)) to be approximated by the 
parabola V — ^rri^clP. While this might significantly affect 
the actual motion obtained for 4>{t), it is irrelevant for our pur- 
pose of assessing the effect of the extra adiabatic approxima- 
tion. 

The results of this numerical comparison of the nonequi- 
librium and adiabatic approximations are broadly similar to 
those of the parallel calculations reported in lll3ll for an in- 
flaton decaying into scalars. In both cases, we find exam- 
ples of conditions for which the adiabatic approximation pre- 
dicts overdamped motion (and so might seem, at least on the 
grounds of self-consistency, to be reliable), while the motion 
resulting from the "full" nonequilibrium evolution is oscilla- 
tory, with relatively little damping. (In neither case did we 
find conditions under which the frictional effect generated by 
the nonequilibrium evolution is sufficient to cause overdamp- 
ing; however, we were unable to explore the large parameter 
space exhaustively.) Even with the adiabatic approximation, 
underdamped motion seems to be more typical of the model 
we studied, but the magnitude of the frictional effect seems 
to be significantly overestimated by this approximation. For 
the model studied in this work, indeed, it seems that friction is 
always overestimated by the adiabatic approximation. In the 
case of inflaton decay into scalars, the adiabatic approxima- 
tion yields a friction coefficient that contains a factor cfP, and 
appears to overestimate the frictional effect when (p is large, 
while underestimating it when (j> is small. 

Beyond the formal difficulty discussed in section UlICI the 



reasons for the failure of the adiabatic approximations are hard 
to identify in any detail. Poor thermalization may be an issue, 
but as noted earlier, it seems implausible that the large dis- 
crepancies we observe are due merely to the quantitative effect 
of number densities that differ somewhat from a true thermal 
distribution. A more plausible possibility is that the discrep- 
ancies could be lessened simply by retaining more terms in 
the time-derivative expansion. Comparison of (|62] | and ( l63b 
makes it rather obvious that, if the adiabatic motion is fric- 
tion limited, so that cj) « —U{(t>)/rj{(f)) and « 0, and if 
the same value of and the same number densities are used 
to evaluate both the adiabatic and the "full" equation of mo- 
tion, a large discrepancy is more or less inevitable. It could 
be, however, that this friction-limited motion does not lead to 
negligible values for higher-order terms in the time-derivative 
expansion, and that this accounts for the discrepancy. Unfor- 
tunately, pursuing the adiabatic approximation to higher or- 
ders is rather complicated (and not entirely unambiguous, in 
view of the evolving temperature we employed); and to our 
knowledge, approximation schemes of this kind are not com- 
monly used in the cosmological literature. 

The possibility of describing the dissipation of inflaton en- 
ergy by a local friction term is primarily of interest in con- 
nection with warm inflation scenarios, and it should be em- 
phasized that the results reported here and in ifTsIl do not in 
themselves determine whether or not such scenarios may be 
viable. To the extent that any generic result can be extracted, 
it seems that the adiabatic approximation overestimates the 
effects of dissipation: in particular, overdamped motion (and 
hence warm inflation) is apparently harder to achieve than the 
adiabatic approximation suggests. However, this may be true 
only of the simple models we have investigated. The some- 
what discouraging conclusion is that no simple analytical ap- 
proximation known to us represents dissipation reliably, and 
numerical investigations of the kind reported here offer a very 
inefficient tool for inflationary model building. 
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